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Abstract The magnetic field structure and the colli- 
mated outflow of rapidly rotating Young Stellar Objects 
(YSOs) are calculated from the stellar source to the 
asymptotic region. The calculations are based on ideal 
MHD and the further simplifying assumptions of station- 
arity and axisymmetry. The star-disk-jet system can be 
completely defined by the Grad-Shafranov (GSS) equa- 
tion, describing the structure of the magnetospheres, and 
the wind equation, which are given by equilibrium per- 
pendicular and parallel to the field lines. Both equations 
must be solved simultaneously to obtain a self-consistent 
solution. General solutions of the Grad-Shafranov equa- 
tion are not yet available. Here we discuss an analytical 
model for the magnetic flux surfaces which is a solution for 
small and large radii. This model assumes a stellar dipole 
field and a gap between the star and a disk at a distance of 
a few stellar radii. Due to the features of the disk no field 
can penetrate the disk and the resulting opening of the 
field lines close to the polar cap is obtained as a compu- 
tational result. In addition, our model guarantees the col- 
limation of the outflow into a cylindrical shape at asymp- 
totic jet radii of several thousand stellar radii. This model 
for collimated outflows reproduces all essential properties 
of magnetospheres for rapidly rotating stars. The result is 
used as input to the wind equation. This problem is com- 
pletely integrable, determined by five constants of motion: 
the total energy E, the total angular momentum L, the to- 
tal mass flux ij along a flux surface, and the total entropy 
S in a flux surface, together with the rotation fl^ of the 
field lines. For adiabatic plasma flows, this problem is al- 
gebraic and can easily be solved. Pressure is neglected in 
our computation. 

The theory of axisymmetric magnetospheres around 
rapidly rotating stellar sources is outlined including elec- 
tric flelds due to the rapid rotation. Gravity of the central 
object is consistently built into this theory. Due to the in- 
jection of plasma either from the stellar surface or by in- 
teraction with a surrounding disk, these magnetospheres 
are neither vacuum solutions, nor force-free. A consistent 
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wind theory is developed which contains the Newtonian 
theory as a classical limit. 

Current-carrying plasma flows will lead to a collima- 
tion of the magnetospheric structure into a cylindrical 
shape. Particular solutions are discussed for the asymp- 
totic collimation. We show that the asymptotic structure 
is essentially determined by electric forces in the pinch 
equation, and not by centrifugal and pressure forces. 

Key words: Magnetohydrodynamics (MHD) - stars: pre- 
main sequence - stars: mass loss - ISM: jets and outflows 



1. Introduction 

The phenomenon of highly collimated high- velocity bipo- 
lar outflows occurs frequently in Astrophysics for a large 
range of objects such as protostars and Active Galac- 
tic Nuclei. The formation of jets seems to be closely re- 
lated to the presence of magnetic fields and the existence 
of a gaseous disk around a rapidly rotating central ob- 
ject (Camenzind 1990). It is widely believed that due to 
these rotating magnetospheres, the mass flow can be ac- 
celerated and collimated on scales of a few tens of AU 
as shown by HST observations of Young Stellar Objects 
(Ray et al. 1996). Nevertheless, their origin is still under 
discussion. Protostellar jets are most probably disk winds 
or stellar winds. The structure of these magnetospheres 
follows from solutions of the Grad-Shafranov equation for 
force-balance perpendicular to the field lines. Magneto- 
spheres generated by isolated stars are typically three di- 
mensional objects (e.g., the Earth's magnetosphere, pulsar 
magnetospheres or stellar magnetospheres) and as such 
are very complicated to model consistently. There are a 
few types of systems where the magnetosphere is con- 
strained to be axisymmetric. Accretion disks will produce 
highly axisymmetric structures and stars in close interac- 
tion with accretion disks are also sources of axisymmetric 
magnetospheres. Magnetized rotators are natural sources 
of energetic plasma fiows. We investigate the dynamics of 
the wind in these collimated magnetospheres for axisym- 
metric and stationary flows. The shape of the magnetic 
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flux function determines whether the wind is accelerated 
or decelerated. 

In this paper we develop a general theory for steady ax- 
isymmetric ideal plasma flows driven by rapidly rotating 
magnetospheres of stellar sources including contributions 
from electric fields. This theory also includes gravity, such 
that outflows driven by neutron stars and black holes are 
covered. With the elements developed in this paper one 
can calculate plasma flows driven by young steUar sources, 
where the central object is rotating near its breakup ve- 
locity, as well as winds driven away by rapidly rotating 
neutron stars. This is an extension of the work by Mo- 
barry & Lovelace (1986) including new insights into this 
problem from the last ten years. 

An overview of the attempts to approach the prob- 
lem of solving the Newtonian Grad-Shafranov equation 



and the wind equa tion can be found in Lery et al. (1998) 
Lcry et al. (1998) also present a model for stationary, ax- 



isymmetric MHD winds without the need for self-similar 
assumptions as is done here. Whereas our model gives 
the structure of the magnetospheres for all regions, they 
find a solution of the Grad-Shafranov equation for the 
Alfven surface from the Alfven regularity condition and 
at the base of the flow. Their magnetic field structure 
for small radii is conical in shape, which does not take 
a stellar dipole field or the influence of a disk into con- 
sideration. We find that the shape of the flux surfaces 
close to the stellar source is important for the accelera- 
tion of the plasma. The Alfven point and the fast mag- 
netosonic point are rather close to each other. The fast 
magnetosonic surface lies at radii smaller than the inner 
disk radius. The coUimation of the magnetospheres into 
a cylindrica l outflow is obtain ed by imposing an external 
pressure in Lery et al. (1999j| . That is not necessary for 
our model. When plasma flows reach radii near the light 
cylinder, electric forces become important for the mod- 
elling of the underlying magnetosphere for fast rotators 
(Camenzind 1986a; Fendt et al. 1995; Fendt & Camen- 
zind 1996). This does not mean that plasma flows are ac- 
celerated to relativistic energies, only that electric fields 
contribute to the forces as much as magnetic fields do. The 
electric field can only be neglected, at radii much smaller 
than the light cylinder radius (at least only a few percent 
of the light cylinder) . The poloidal current follows from the 
MHD wind theory. For [Lery et al. (1999)| the li ght cylin- 
der is of no importance. In Lery et al. (1998) the wind 
equation is solved from the source to the fast magnetosonic 
point, whereas our wind solution is global and goes up to 
the asymptotic region of the jet. In contrary to our cold 
wind solution they use a polytropic equation of state. For 
a cold wind our fast magnetosonic surface does not go to 
infinity, because we impose a dipole-like structure on the 
flux surfaces instead of a monopole-like one. In addition, 
their wind solution represents the Newtonian limit of our 
derivation. Newtonian MHD is perfect for modelling the 
magnetic structure of the solar wind (Mac Gregor 1996). 



But winds ejected from protostars rapidly rotating with 
periods of a few days could pass the light cylinder given 
as (Camenzind 1997) 

(1) 

Such a radius is interestingly close to the observed coUi- 
mation radii. While the light cylinder is not essential for 
the plasma kinematics, it does influence the collimation 
processes at large radii by means of electric fields due to 
field line rotation. 

The paper is structured as follows. In Sect. |2| a consi- 
stent wind theory for rapidly rotating objects is presented. 
Axisymmetric fields and stationary flows are investigated. 
The Newtonian case follows quite naturally as a limit of 
the special relativistic MHD. The polytropic equation of 
state is included. A general expression for the wind equa- 
tion is given and the special solution for a cold wind is 
discussed. In Sect. ^ the Grad-Shafranov equation is stud- 
ied. Sect. ^ gives a discussion of the collimation processes 
for spherical outflows. The importance of the electric flelds 
close to the light cylinder is shown. In Sect. || a model is 
presented to reproduce the characteristics of collimated 
outflows near to and far from the star. A calculation of 
the magnetic flux surfaces and the resulting wind is given 
in Sect. ^ for a typical young, low-mass star. Also the in- 
fluence of various parameters of the model is studied. We 
draw our conclusions in the last section. 



2. A consistent wind theory for rapidly rotating 
stars 

A modern treatment of magnetohydrodynamics includ- 
ing the effects of gravity must be based on a general 
relativistic approach in the 3-1-1 split of space-time (see 
e.g. Thorne et al. 1986; Beskin & Pariev 1993; Camenzind 



1996). When the rotation of the central object is not im- 
portant, the line element is given by the Schwarzschild 
decomposition 



ds^ = -a^ ^2 ^^2 ^ J_ ^^2 ^ ^2 ^^2 ^2 gj^2 g ^2 ^ 



Here a = a/1 — 2GM*/c^r is the redshift factor which 
accounts for the gravitational redshift between a local ob- 
server and infinity. It describes the effects of gravity due 
to the presence of the central star of mass Af*. In the 
Newtonian limit we have 



ac^l + $/c^ 



(3) 



where $ is the Newtonian potential. In the following we 
also use the cylindrical radius R — rsinO. The time t is 
the global time. 

The spin J, of the central object is also a source of 
gravity in general relativity (the so-called gravitomagnetic 
effect). This spin would produce an off-diagonal element 
gtif, in the above line element, g^^ oc J^,/r^, which can 
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be important for neutron stars rotating with millisecond 
periods (Camenzind 1986b). This effect is crucial for the 
treatment of plasma flows on rapidly rotating black holes 
(Camenzind 1996). Since we are interested here for appli- 
cations to stellar sources, we will neglect these effects. 

Local fiducial observers (FIDOs) are represented by 
orthonormal tetrads 



O 
eg 



1 

ac 
r 



dt 



= adr , 



(4) 
(5) 



2.1. 3+1 split of electrodynamics 



On a curved spacetime background electromagnetic fields 
must be given in terms of the Faraday tensor F'^'^ , and 
Maxwell's equations are given in tensorial form. With re- 
spect to the particular observer field defined above, we 
can however split the Faraday tensor as in special relativ- 
ity into an electric and magnetic part 



Similarly, we get a charge density 

Pe= O J 

and a current density]^ given by 
J ^ PeO + J . 



(6) 



(7) 



(8) 



In terms of the electric fields E and magnetic 
fields B, Maxwell's equations assume the familiar form 



(Thorne et al. 1986) 



47r 



V • £; = 47rpe , 

V B = 0, 

Charge conservation is given as 
0. 



aj. 



(9) 
(10) 

(11) 

(12) 



(13) 



In this treatment, the redshift factor a only appears 
as a form factor, but there are no additional couplings be- 
tween electromagnetic fields and gravitational forces. This 
is in contrast to the treatment on rapidly rotating back- 
grounds, where the gravitomagnetic force couples directly 
into Maxwell's equations (Camenzind 1996). 



^ We use the notation of Thorne et al. 1986 for the four- 
vectors. 



2.2. Axisymmetric fields 

As in the Newtonian theory of plasma confinement, the 
above equations can be considerably simplified for axisym- 
metric structures. Here, we can define the poloidal flux 
function ^'(r, 9) 



"^{r.e)^^ [ Bp-dA, 

J A 



(14) 



and the axial current flowing downward along the central 



/(r, 



ajp ■ dA = — - uRBt ■ 



(15) 



A denotes an upward directed surface around the central 
axis. The level surfaces of are the magnetic surfaces 
which are generated by rotating field lines. The level sur- 
faces of I are the surfaces where the poloidal currents flow. 

Due to stationarity and axisymmetry, the electromag- 
netic fields assume a simple form 



B = 4 V^- X 

-ft 



21 

Rac 



Er. 



Et^O. 



ac 



(16) 
(17) 



Infinite conductivity requires B ■ E = 0, and thus the in- 
tegration constant ri^(5') represents the angular velocity 
of field lines. This is not the angular velocity of a plasma 
particle. 

2.3. Plasma motion 

We t reat plasma m otion only in the one-fluid approach 
(see Khanna 1998 ). For a one-component non-viscous 
plasma the energy-momentum tensor has the form 



T={p + P)u(g)u + Pg + T^, 



(18) 



p is the total energy of the plasma including internal 
energy, P the proper pressure, u the four-velocity of 
the plasma and g the metric tensor. We decompose the 
energy-momentum tensor T into horizontal and vertical 
components with respect to fiducial observers 



T = eO<^0 + 0(» S+ t . 



(19) 



e is the energy-density with respect to local observers, S 
the momentum flux and t the stress tensor. 

With the decomposition of u = 7(0 + 1;) we obtain the 
following components with respect to FIDOs 



e = 7^ (p + Pv^ 



(20) 
(21) 



S = (p + P) + Sen-, = POMT^ V + Sam, 

T = {p + P)-f'^V(^V + Pli + Tom 

= S^-u + Py+Tom • (22) 
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Here, po is the rest mass density, p, — {p + P)/n the rel- 
ativistic specific enthalpy and 7 the Lorentz factor mea- 
sured by fiducial observers. Sem represents the Poynting 

fiux measured by FIDOs and t cm the Maxwell stresses 

1 



— LE" 

Stt V 

— ExB. 

An 

1 , 

An 



E~ Bi 



B + -g (E^ + B^) 



(23) 
(24) 
(25) 



Using the 3+1 split of the affine connection, one can now 
derive the 3+1-split of the hydrodynamic equations V • 
T = 0. Energy conservation, given by O • (V • T) = 0, can 
be written as (Thorne et al. 1986) 



1 



dte 



1 



V-{a^S). 



(26) 



de 

dr a a' 

The right-hand side in the energy conservation is the fa- 
miliar divergence of the energy flux, with one factor of a 
inside the divergence to account for the gravitational red- 
shift of the energy and the other to convert the differential 
proper time element in the definition of the flux into to a 
differential universal time element. 

Similarly, Euler's equations, given by /i • (V • T) = 0, 
where h is the projector into the plasma rest frame, assume 
the form 



7 o -1 -1 

— = -dtS= -eV(lna) - - V • (a t). 
dr a a 



(27) 



The quantity —Vina represents the local gravitational 
force measured by observers. 

2.4- Stationary flows 

The plasma flow is given by the continuity equation 

V • (ajnv) = . (28) 

Similarly to the Newtonian theory, one also obtains for 
the poloidal velocity 



B 



p ' 



(29) 



and, from the frozen-in condition, that the plasma is 
guided by the structure of the axisymmetric fiux surfaces 
5* = const 





60 . 


an 


a 



(30) 



As a consequence of Faraday's induction law, the rotation 
of the field lines remains constant at il^(^'). The above 
relation tells us that the particle fiux in a flux surface 
remains constant. 

For stationary and axisymmetric flows the energy 
equation ( p6|) is now simply given as 

V • (a^S) = 0, 



and similarly for the momentum equation (|2 

-V • (a T) + eVlna = 0. (32) 
a 

By means of the MHD relations the energy flux can be 
written as 

S^fij^Bp + —-^Bp. (33) 

In the stationary case, this provides the following expres- 
sion for the energy equation 



V • fijarj Bp 



■ -Bp I = . 



(34) 



Using fiux conservation V • Bp = and {Bp ■ V)!] = we 
find conservation for the total energy per particle 



£;(*) = p.aj + 



2TTr] 



(35) 



A similar conservation law follows from the angular 
momentum equation 

e0-V-(aT)=O. (36) 
Due to axisymmetry, all other terms drop out. Since 

60- t — tpj, — fiUpj — — BpBff, 
47r 



^Bpj + —Bp 



wc get 



V • i^miBp + ^Bpj = , 

or the total angular momentum conservation 



L(*) = ^,J 



I 

27r77 ' 



(37) 



(38) 



(39) 



with the specific angular momentum j = R^fl/a^. In the 
total energy, the kinetic part is redshifted. Together with 
the above relation (^^ for the toroidal motion, 

u, = ^B, + ^, (40) 
an a 

this provides three equations for the unknowns 07, j (or 
u^) and /. 

In addition to these relations, an adequate equation of 
state is needed. For n = n{P, s) and T ~ T{P, s) the first 
law of thermodynamics implies 

d^i=-dP + Tds. (41) 

In particular, for P — KQ{s)n^ with a polytropic index 
r, one finds for the specific enthalpy 



(31) f^ = mc + 



T - 1 



Ko{s) n 



r-i 



(42) 
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With this expression, we can derive the Newtonian hmit 
for the total energy, a — 1 + <&/c^, 7=1 + 



1 



r p 



(43) 



E = mC H — V . , ^ , 

2 r - 1 po 27r?7 

which is the standard energy conservation used in Newto- 
nian MHD. 



2.5. Lorentz factor, angular momentum and poloidal 
current 

For given integrals of motion L{'^), il^(\E'), ?/(^') 

and i4r(\[') the conservation laws for stationary and ax- 
isymmetric plasma flows can be reduced to give relations 
for the Lorentz factor, the angular momentum and the 
poloidal current flux function I{R, 5") 



= Q!7(\l/,i?) 



E a^(l-e)-M^ 



(44) 



/(*,i?) 



2Tn]E a^e-{R/RLf 



nP - {R/RlY - M'^ ' 
where the quantity M, defined as 



NP 



v2 ,,2 



ip Ct Up X 



(46) 



(47) 



represents the (poloidal) Alfven Mach number. We nor- 
malize radii in units of the light cylinder R^ , x = R/Rl- 
The redshift factor a corrects for the singular behaviour 
of the poloidal velocity at the horizon of the black hole. Q 
The Alfven Mach number M is a well-behaved quantity, 
e.g. at the horizon. The quantity e defined by 



e(*) 



E 



(48) 



is a measure for the amount of energy carried by the 
electromagnetic fields. In special relativity, e < 1, and 
e = R\/R\ is a direct measure of the position of the 
Alfven surface in terms of the light cylinder radius Rl. 
For relativistic flows, Ra Rl, and therefore e — > 1. For 
the solar wind ~ 10~* and for winds ejected by rapidly 
rotating young stellar objects we find ~ 10~^. In the 
case of the Crab pulsar ~ 0.999 and the wind dynamics 
is essentially relativistic. 

This quantity e determines the position of the Alfven 
surface (positions, where the denominators in the above 
equations vanish) 



a\RA)~{RA/RLf ^Ml. 



(49) 



This equation has in general, for given rotation fi^, only 
one solution, the outer Alfven surface corresponding to 

^ These relations have been derived for the first time in Ca- 
menzind (1986b) without using 3-1-1 split. 



the special relativistic one. The inner Alfven surface only 
exists in black hole magnetospheres. The light cylinder 
surfaces are special solutions of this equation for M\ = 0. 

The regularity condition at the Alfven point also de- 
termines the total angular momentum L as a function of 
the position of the Alfven point 



M 



#2 



L 

E 



RW (1 - e) . 



(50) 



which becomes in the Newtonian case simply L = R\Q^ . 
2.6. The Newtonian limit 

There is much confusion in the literature as to the Newto- 
nian limit of the special relativistic MHD. The two equa- 
tions (^) and (|4^) are well-known in Newtonian MHD 
(a = 1, e <C 1 and Rl — > 00), where they form the ba- 
sis for a treatment of axisymmetric MHD winds in stellar 



systems (Beyvaerts & Norman 1996) 



RB^ 



R^VL 



L 



R2^F _ ^2^ 

1-M2 

R^n^ 



l-AP 



(51) 
(52) 



In Newtonian MHD, the Mach number is simply AP = 
Vp/v\. The first equation indicates that plasma is corotat- 
ing with the field lines for low Mach numbers, ^ 1, 
and that the specific angular momentum j is equal to 
the total angular momentum L for high Mach numbers, 
^ 1. These equations tell us that Newtonian MHD is 
only valid inside the light cylinder Rl — c/ft^ . Since for 
stellar objects, e <C 1, a comparison between Eqs. (^5|-^6|) 
and (^l|-|5^) tells us that the Newtonian limit is correct 
for the wind equation, as long as ^ x^ for all relevant 
radii. This approximation is therefore not justified e.g. for 
rapidly rotating magnetic surfaces generated by accretion 
disks around black holes. Even in the case of magnetic 
surfaces generated by rapidly rotating young stars slight 
modifications occur in the wind equation. 

In the Newtonian limit, the energy conservation at the 
injection point can be written as 



E 



- 1 = 



2 c2 



r 



r — 1 mc?7i.f 



(53) 



i'Kfnrjc^ 

Using the solution for the current at the injection point, 
xl < e, < 1 and = kT^/mP, 

2'Kfiri e 



I 



1 



(54) 



this shows the various contributions to the total energy 
-B _ 1 V? r e 



" c2 + 2 c2 



r - 1 



p* 



1 



(55) 



The parameter e is therefore a real measure of the Poynt- 
ing flux at the injection point in units of rest mass en- 
ergy fiux. For the solar wind one finds e© = 10~^ and 
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= 10~^, and therefore the solar wind is thermally 
driven and not Poynting flux driven (see also MacGregor 
1996). For rapidly rotating protostars, the Poynting flux 
increases, e ~ 10^^ which demonstrates that these winds 
are Poynting flux driven, even if they would be ejected 
from a hot corona with ~ 10^^. 

The ratio e/p, is the true measure for the question 
of Poynting flux driven outflows. For relativistic winds, 
e ^ 1, and the Poynting flux can then even exceed the rest 
mass energy flux. 

2. 7. The wind equation 

In ideal MHD, plasma flows along the magnetic surfaces 
5* = const. Since the Lorentz factor, the angular momen- 
tum (or angular velocity 11) and the poloidal current func- 
tion / are essentially only functions of the radius along the 
flux surface and the Mach number, the normalisation of 
the plasma 4-velocity, UaU°' = — 1, 



7^ — = 1 , U — "fV 

leads to the equation 



1 2 J 



(56) 



(57) 



Using the solutions for 7 and j, this can be rearranged 
into the form 



(58) 



where the subscript S is for 'Schwarzschild'. With radii 
measured in units of the light cylinder radii — c/il^ , 
the potential function Vs can be expressed as 



Vsix;M^,e) 



with the definitions 

Ds = a^- {R/RlY - A/2 , 

Fs = x^[a^{l-e)-M^Y 
- a\l-e)x^-M^e]\ 



(59) 



(60) 



(61) 



This relation is also correct near the Schwarzschild hori- 
zon, since aup remains finite there. 

The above equation represents an algebraic relation 
for for each position along a flux surface. This follows 
from the relation between u„ and Br, 



(aupf = 



o2 

^ M\ 



Multiplying the equation by x"* we obtain 



2 Fsix;M',e) ( E\ ^ . 



Bix'^ 



(62) 



(63) 



The last term can be written in terms of the flux function 
BnR^ 



<i>^ix) = 



Bp^R'^ 



(64) 



which is normalized by the footpoint i?* and its magnetic 
field strength Bp^. In addition, we introduce the dimen- 
sionless magnetisation parameter 



{Bp,Rl){^)c 
^^:^l7^{^)Rl{■^) 



(65) 



which is the analogue of the magnetisation parameter in- 
troduced in the Newtonian discussion. This altogether 
leads to the equation 



Fs{x-M\e) (E 



Dl{x-M^) \^i 



-) ^a^x' + al'^lM^ 



(66) 



This fundamental equation clearly shows that the solu- 
tions of the wind equation, when formulated for the Mach 
number M, depends on the following parameters: 

— the parameter e, which defines the position of the 
Alfven point 

— the dimensionless energy E = E /mc? ^ 

— the magnetisation parameter cr*, 

— the flux tube function function $iif 

— two additional parameters hidden in the specific en- 
thalpy 



mc 



F - 1 



Ml 
M2 



r-i' 



(67) 



These are essentially 5 parameters for each flux sur- 
face ^' = const. It can be shown that the relativistic 
wind equation has the same critical points as the non- 
relativistic one: the slow magnetosonic point, the Alfven 
point Ds{Ra) = and the fast magnetosonic point. It is 
very important in this regard that the light cylinder is not 
a critical point |^ 

The requirement that a wind solution passes through 
all three critical points fixes therefore three of the five pa- 
rameters. We may consider cr, and p^ as free parameters 
which are fixed by injection physics. 

When pressure is as important as Poynting flux, p^ ~ 
e, we can apply the expansion of the wind equation 



x'^FsE^ = [D^a'^x^ + al'^^D^M'^) 
with p. ~ ^/mc^ 



2F 



F- VM2 



Mi 



0{pl). 



(68) 



(69) 



These points are modified somewhat when the structure 
of the magnetosphere is taken into account (Tsinganos et al. 
1996). 
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Since protostellar outflows are probably nearly isothermal 
(Paatz & Camenzind 1996), F ~ 1.2, the energy equation 
can also be transformed into a polynomial form for the 
square of the Mach number 



m2„M 



2n 



0, 



(70) 



with a suitable power Up {up = 9 for T — 1.2). This form 
of the wind equation can then be used to study all real 
branches (for a Newtonian study, see Paatz & Camenzind 
1996). 



2.7.1. The cold wind equation 

In the cold limit, <C e, we may neglect the pressure 
terms in the specific enthalpy, /i ~ rnc^, and the wind 
equation then provides a 4th degree polynomial in the 
square of the Mach number for the poloidal velocity 



(71) 



with the following polynomial coefficients 

mo = E-'a'il - efx^ia" - x") - a^x^a^ - x^ f , (72) 

ma = ~2E^a^{l~ e^x'^ + 2a^x^{a^ - x^) , (73) 

m4 = E^x^ix^ ~ o?^) - c?x^ - al<^''{o? - x^ f , (74) 

me = 2f72$2(c,2 _ 2.2) ^ (75) 

(76) 



and with ra2n — "t,2„(£', e, cr,, $; x). Without gravity 
(a — 1), these expressions are identical with the flat space 
expansion (Camenzind 1986b). In the case of a cold wind 
'Wp(xinj) = 0, where Xinj is the injection point of the plasma 
from the star, and mo(a:inj) = 0. Consequently, we get a 
relation between injection radius, energy and Alfven ra- 
dius 



x\^\- 



Ea 



(77) 



hence m2n = m2„ (xinj , cr* , $; x) . A cold wind has only two 
critical points, the fast magnetosonic one xfm, where the 
wind velocity must equal the fast magnetosonic velocity, 
and the Alfven point. The slow magnetosonic point xsm 
disappears, with the consequence that the mass flux or the 
magnetisation along a flux surface becomes a free param- 
eter. The solution of Eq. (^) is uniquely determined by 
the critical energy £'crit(^)- Eq- ( fn] ) is a 4th order poly- 
nomial. For this reason the wind equation has 4 solutions 
Up,i{x), {i = 1, ...,4) at each radius x. We consider nega- 
tive or complex solutions to be 'unphysical'. We obtain two 
positive branches of which the 'physical' solution is the one 
which passes through the critical points and starts with 
zero poloidal velocity. In order to determine the critical 
solution, we use a bracketing method varying the energy 



parameter E{^^) until the critical solution is found. For 
energies E < E'crit, the wind ceases before being able to 
reach the fast magnetosonic point or the asymptotic ra- 
dius, hence there is no continuous wind solution for all x. 
In the case E > Edit , the two branches of the solution do 
not meet at the fast magnetosonic point (see also Fendt 
& Camenzind 1996) and no wind solution is found cither. 
We fix the free parameter e by choosing the values for the 
injection radius a;inj and the magnetisation parameter cr*. 
For small a.^, the energy must be determined up to the 
order 10^^ (t^. 

3. The transfield equation and its action integral 

Ampere's equation determines the form of the magnetic 
fiux surfaces, known as the Grad-Schliiter-Shafranov 
equation for ^! 



a \B? 



A*^ EE — V • f -^V^- 



47r 



(78) 



V is here the covariant derivative with respect to 3-space. 
This divergence-operator can explicitly be written as 



~ drV" dr)^ r2 86 I sinff 39 



(79) 



and is equivalent to the fiat space analog except for the 
redshift factor a. 

The split-monopole 

= l-^l - cose*) (80) 
is still a solution of the vacuum equation. A dipole 



sm 



(81) 



would be modified somewhat by the gravitational poten- 
tial. 



3.1. The Grad-Schliiter-Shafranov equation 

When the magnetosphere is filled with plasma, currents 
will deform the vacuum solutions. The current density 
follows from the transverse component of the equation of 
motion 



eVj.a + V ■ {a t)j_ = 0. 



(82) 



The explicit form of this equation has been derived in var- 
ious papers (Mobarry & Lovelace 1986; Camenzind 1987; 
Beskin et al. 1995; Fendt et al. 1995) 

^'V- (^V*) =5(/,m2,...) (83) 

with Ds as defined in Eq. ( |60| ) and g{I,M'^, ...) a source 
function. 

This equation has various interesting limits. One is 
the force-free limit with vanishing pressure and Af ^ = 



8 



E. Broitmoser & M. Camenzind: Collimated outflows of young stellar objects 



(Blandford & Znajek 1976; Mobarry 
Okamoto 1992; Fendt 1997) 



Lovelace 1986; 4. Collimation of spherical outflows 



Ds 



a 



= 0. 



(84) 



In this case, Ds — c? — and the source of the mag- 
netic flux is only given by the poloidal current Even 
in this limit, the equation is highly non-linear and the 
Alfvcn surface now coincides with the light cylinder. Such 
solutions are interesting for black hole driven outflows. 

Another limit which has been discussed quite often 
in the literature is the Newtonian limit without gravity, 
a = 1 and R <C Rl (Heinemann & Olbert 1978; Pelletier 
& Pudritz 1992; Rosso & Pelletier 1994) 



1 



i?2 



V^- =.9Ar(/,M^...). 



(85) 



Rosso & Pelletier (1994) have used this limit for the dis- 
cussion of disk magnetospheres. They have derived a La- 
grangian description of the form 



(86) 



depending on a potential function 1}q{M'^ ."^^R^z). It is 
then evident that this equation can be derived from the 
Lagrangian 



Lo(M2,*;i?,z) 



1 M2 - 1 



(V*)2 + J7o(Af2,^';i?,z)(87) 



2 i?2 

and the proper action of the system is simply 
^i^-] = J 2TTLo{M^,'b;R,z)RdRdz . 



The above formulation is indeed only a special case of the 
action discussed by Mobarry & Lovelace (1986). 

3.2. On the causal structure of self-consistent 
magnetospheres 

As shown above, the concept of light cylinders appears 
naturally in ideal relativistic MHD. This condition is 
largely satisfied in thin outflowing plasmas, in partic- 
ular when plasma will be reheated by some processes 
(shocks e.g.). The light cylinder is not a critical surface 
of the Grad-Schliiter-Shafranov equation (see e.g. Beskin 
& Pariev 1993) except in the force-free limit. As in the 
Newtonian case, the relativistic Grad-Schliiter-Shafranov 
equation has the same critical surfaces, the slow and the 
fast magnetosonic surface, as well as the Alfven surface 
(Heinemann & Olbert 1978, Tsinganos et al. 1996). The 
system becomes hyperbolic beyond the fast magnetosonic 
surface, so that no back-reaction with the star is allowed 
from this region. 



The collimation of magnetized winds into a conical outflow 
is very similar to the plasma confinement in a z-pinch. In 
this region we can neglect gravity. The equilibrium con- 
dition for a non-rotating axisymmetric current carrying 
pinch configuration (Freidberg 1987) 



47r 



87r 



- Sir 



47ri?3 



(-V^i?) (89) 



has to be generalized to include various relativistic effects 
(Appl & Camenzind 1993a, b) 



Bl 



K [1 - - x"") = (1 - a;2) Vj ^ + V I -r^ + V I P 



47r 



^8vr 
i?3 ~ 47ri?3 



87r 



(-V^i?). (90) 



Pressure gradients acting perpendicular to the magnetic 
surfaces are in equilibrium with electric forces, centrifu- 
gal forces, pinch forces produced by poloidal currents, as 
well as with curvature forces expressed at the left hand 
side. This latter force is modified by the Mach number M 
and the light cylinder, and /i = {p + P) /n. A comparison 
between the Newtonian expansion ( p9| ) and the special 
relativistic expression (|9^) clearly demonstrates the influ- 
ence of the light cylinder in the force-balance equation. 
The major distinction is found in the first expression on 
the right-hand side, (1 — x^) V ±Bp/8Tr. This term changes 
its sign at the position of the light cylinder, a; = 1 , as a 
consequence of the existence of electric fields in Maxwell's 

stress tensor t em , Eq- ( p5|) . Since the observed collimation 
radii are larger than the light cylinder, field lines anchored 
in the star will cross the light cylinder, at least when the 
stellar rotation periods are of the order of a few days. 

The specific angular momentum of the plasma j = 
R^rt = Ru^l^ is now given by the wind theory 



J = Rlc 



(1 - e)x^ - APe 



(91) 



and depends on the Mach number M . Outside the Alfven 
surface, > 1, this can be expressed as 



(1 - e)x^ 
AP 



(92) 



Outside the Alfven surface, the specific angular momen- 
tum is essentially constant and has the value 



,2 nF 



j ~ €Rlc = R\n 



(93) 



Similarly, the pinch current / is not a free quantity, 
but also follows from the wind theory 



2'Kr]E e — x"^ 

1 - a;2 - M2 



(94) 



E. Broitmoser & M. Camenzind: Collimated outflows of young steUar objects 



9 



Using the expression for the Lorentz factor, the total en- 
ergy E can be replaced by, a = 1, 



27r r//i 

r 



M2 



(95) 



For radii > e and Mach numbers Af^ > 1 , i.e. es- 
sentially outside the Alfven surface, the current can be 
approximated by 



(96) 



With the definition of 77 this gives the universal expression 
for the current 



5 S 



(97) 



always valid outside the Alfven surface. f3p = Vp/c is the 
dimensionless poloidal velocity of the plasma. 

4-1- Asymptotic collimation 

We discuss in the following the solutions for a cylindri- 
cal pinch, where the curvature forces vanish, k = 0. In 
this case, the equilibrium condition reduces to its one- 
dimensional form 

dR 8tt 



(1 



dR 8tt 



dP 
dR ~ 
P 



2ttRI 



= 0. 



(98) 



The pressure of the toroidal field and the pinch force can 
be combined into one single expression. In addition, we 
normalize the poloidal magnetic field in units of its value 
on the central axis, Bq, y = Bp/B^, 



[l-x') 



dy_ 

dx 



-4:xy + 



1 dP Stt dP STTfin f 



dx 



Bl dx 



RlBl x3 



0.(99) 



Electric fields, Eji_ — {R/ Rl)Bp, due to the rotation of the 
magnetic surfaces modify the classical pinch equilibrium 
in two ways. First, the equation has a singular point at 
the light cylinder, and secondly, the additional term 4xy 
is crucial for the form of the equilibrium. In contrast to 
force-free models (Appl & Camenzind 1993b), the specific 
angular momentum j of the plasma and the poloidal cur- 
rent / are not arbitrary, but follow from the MHD wind 
theory. 

In particular, the last term can be written as 



,100) 

Rf^B^ x-^ Sq Rf x-^ Rf x'^ 



where Va is the local Alfven speed and v^^a = Ra^^ the 
rotational velocity at the Alfven point. We also expressed 
the specific angular momentum in terms of its character- 
istic value, j — j R\il^ . This term contributes therefore 



to the second term in the pinch equation (|99|). This shows 
therefore that the centrifugal force is not important near 
the light cylinder and beyond, whenever the Alfven radius 
is far inside the light cylinder. In this domain the electric 
fields provide the important forces. Far inside the light 
cylinder, however, centrifugal forces are important. 

4-2. Pinch solutions without angular momentum 

We first discuss solutions of the relativistic pinch equilib- 
rium neglecting the influence of the centrifugal force. In 
this case the equilibrium is determined by the condition 

(1 2) 4 ^_ 1 + - 

dx x^ dx Bq dx 

The toroidal current follows to a good approximation 



(101) 



I ~-n^{R'^Bp)^ 



to the equation, 7^ = 1/ — /S^ 
^ dy , 4 



1 



-cBoRl (x'^) — 



(102) 



1 

P~p — - 

and we prescribe the plasma pressure in terms of P = 
PBp/SiT with a constant plasma beta. This leads finally 



1 + /3 + 



dx 



xy — . 



(103) 



ip^p / ^pPp 

This homogeneous equation has the remarkably simple so- 
lution 

C 

yi^) = n^R^ 2/2^2 (104) 

[1 + p + x'^/x^y 

with 

xc = jp Pp ^ 10-^ , Rc ^ jpl3p Rl ^ 10'^ Rl (105) 

as the core radius of the jet with l3p ~ 10^^ and 7^ = 1 
for protostellar outflows. Inside the core, given by the light 
cylinder radius and the poloidal velocity /3p, the poloidal 
magnetic field is essentially homogeneous. It decays as a 
monopole field outside the core. Bp cx 1/i?^ for R^ R^. 

This is not unexpected since the jet is formed by a su- 
permagnetosonic wind. The current function / increases 
quadratically inside the core and saturates towards a con- 
stant value outside the core, loc — cBqRlPp- This is ex- 
actly the structure of the force-free solutions found previ- 
ously (Appl & Camenzind 1993b). 

The flux function has the same form as in the force- 
free solutions 

* cx ln[l -I- /3 + (i?/i?c)^] ■ (107) 

This is not unexpected, since in the collimated regime the 
magnetic structure is force-free, / = 1(5'). 

MHD jets always have a core-envelope structure, along 
the axis the magnetic field is predominantly longitudinal, 
beyond the core it becomes dominated by the toroidal com- 
ponent. The light cylinder R^ — 100 AU is the basic scale 
for the jets, the core is extremely small, Rc — lOi?* . 
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4-. 3. Pinch solutions including angular momentum 



5.1. The model 



In this section we investigate the behaviour of the pinch 
equation inside the Ught cyhnder including centrifugal 
forces 



dx 



y 



1 dl^ Stt dP _ ^ 
dx Bn dx 



(108) 



The centrifugal force must be regularized inside the Alfven 
surface. A suitable function which guarantees corotation 
inside the Alfven point is given by 



[R/Ra? 
1 + [R/RaY 



(109) 



Using the expression for the current and neglecting the 
pressure 



1 



x"^ \ dy 



4 + 1 y = 0, (110) 



^IPI) dx Wpf3l 
with solutions as obtained by a simple integration 



In 2/ 



{Al-fll3l)x + CV{xyp/x^ 



dx . 



(Ill) 



or 



In y = In 



1 + xyxi 



^■^[l + x^/x-^ 



■dx. (112) 



The second integral only contributes in the region of the 
Alfven radius. For the case Rc = Ra and V = const, the 
solution is simply given by 



In y = In 



1 



C Rl 



1 



l + x-^/xl) 4 R\ (l + i?7i?i)2 



-Ci.(113) 



Since C oc R\/ R\, angular momentum modifies the solu- 
tion only near the Alfven surface. These effects drop out 
very rapidly beyond the Alfven radius, and we are left 
with the solution found previously. 



5. A model for collimated outflows 

Fully self-consistent solutions of the Grad-Shafranov 
equation are not yet available. We want to demonstrate 
the main features of collimated outflows in terms of a 
model which satisfies the inner boundary conditions and 
the asymptotic flux distribution for cylindrical outflows. 
The existence of a disk surrounding the star will mod- 
ify the dipolar field in the innermost few stellar radii (see 
Fendt et al. 1995). Beyond the inner edge of the disk the 
resulting magnetosphere is however largely spherical and 
cylindrical in the asymptotic region. 



The flux distribution should be a solution of the Grad- 
Shafranov equation. As a flrst approximation we use the 
following ansatz for the flux function, where ^ consists of 
a homogeneous and a special solution 



*(r, 6) = *''°'°(r, e) -f «'^P'='=(r, 6*), 
with 

The vacuum solution (see Appendix A) is given by 



(114) 



(115) 



*'^°'"(p,( 

where 

7 



2/i, p f ip 



cos^e 



+ 7(1-V^)- 2 

5— -arctan- + — , (116) 

P \ 1 ^JJ 



V2\ 



^p^-l + v/(p2 - 1)2 + 4p2 cos2 e 



(117) 



The quantity p is dimensionless due to the normalization 
of the spherical radius r in units of the inner radius of 
the disk Odisk, P '■= ''/adisk- We assume the existence of a 
gap between the star and the disk due to the magnetic 
pressure. The dipole moment is given by /i* = R^B^. 
For our calculations the stellar radius is i?, = ZRq and 
the maximum magnetic field strength on the stellar sur- 
face B^ = 1000 G. The last term in denotes the 
dipole magnetic flux generated by the star itself, while 
the remaining terms denote the flux contribution from the 
screening of the disk (|Kundt fc Robnik 1980| ); [Aly 1980 ; 
paatz fc Camenzind 1996 ). The integration constant -0 



equals the flux on the open fleld lines for the homoge- 
neous solution. For an integration constant ■0 = all flux 
is forced to go through the h ole between star and di sk and 
cannot escape to infinity ( Kundt fc Robnik 198C ). The 
term ■0/^ does not contribute to the magnetic fleld, but 
guarantees \1/ = on the z-axis. In th e collimation region 
we find for large radii from Eq. (116) 



cosO 1 , 

7 = + '^(— )' 

P P* 

*'^°'"(p,0) = ^^0(l-cos0) + O(^). 

TI" fldisk P 



(118) 
(119) 



This means that vl^hom j^gj^aves like a monopolar type func- 
tion for large radii. 

^spcc describes the magnetic flux produced by a 
toroidal current 



^-«P°=(p,0) = A In ( 1 + —(1 - cos( 



i?2 



• (120) 
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Apart from the normalisation factor A, which will be dis- 



cussed in Sect. 3.3, it contains two essential radii. Rr is 



the core radius discussed in the last section, since for large 
radii we find then 



5» 



spec 



In 



1 



sin^ I 



m 



(121) 



as required by cylindrical coUimation. For small radii, the 
cylindrical radius drops out of the expansion and we get 
a monopolar type solution 



In 



^1 <' 



COS( 



(122) 



The parameter a is a kind of jet radius, since the last 
magnetic surface that can be closed with the star is given 
by ln(a/i?f,)^. In the case of Eq. (120), the flux surfaces 
can be parameterized by means of 



(7-1-02(1- cos 6')/i?2 



sm ( 



(123) 



The angle extends from 9 = Q (asymptotic region) to a 
maximal value 0^=,, where the flux surface intersects with 
the stellar surface, r{9(j) — i?, . The flux surface label C 
is therefore given by this angle 



C=l + ||sin2^^ + -^(l-cos0^). 



(124) 



For the total flux function the coUimation part will fi- 
nally always dominate the homogeneous part. Comparing 
Eq. (119) with Eq. (121) we see that the flrst does not de- 
pend on the radius, hence the latter containing the square 
of the radius will determine the solution for sufRciently 
large radii. 



Note furthermore that we need all parts of Eq. (114) 
to reproduce the dipole influence of the star at small radii 
and the coUimation of the field lines for the asymptotic 
radii. Excluding the monopolar type part of Eq. ( |120| ) we 
cannot obtain realistic jet radii and the magnetosphere 
will coUimate much too quickly. v]/sp'='= alone cannot repro- 
duce the influence of the dipole near the star. The surfaces 



of constant flux for the total flux function Eq. (114) are 
calculated with the Newton-Raphson method. The result- 
ing magnetic fields can be found in Appendix B. 

5. 2. The amount of flux through the hole of the disk 

The magnetic flux ^'hoic through the hole of the disk at 
z = consists of the flux belonging to the disk held 
^disk vpspec between star and disk, the whole flux of 
the dipole minus the part of the dipole on the disk 



(Riffcrt 1980) 



5- 



hole 



2a«*p f ,^ cos^ei 

-f{l-1p) 



sm 



■ arctan ■ 



1 



-, 1 



Jp=0 



*Dip - 



7 

^* sm 



p=0 



(125) 



p=i 



For V' = all other terms except for t he flux of th e 
dipole ^pDip and '^''P'"' wiU vanish (see also [Riffcrt 1980|) . 
Then plasma can only escape to infinity due to th e spe cial 
solution. This characteristic is discussed in Sects. 6.3 and 
|6.4| . For tp and A = the amount of flux through 
the hole is determined by the integration constant tp and 
the disk radius Cdisk via ^'hoie oc — V'/«disk- This means 
that for ■0 = and yl = all flux goes through the gap 
between the star and the disk. This amount will decrease 
with increasing ^ and A. 

6. Results and discussion 

We now present solutions for the magnetospheres and the 
resulting winds which are calculated along the coUimating 
jet flux surfaces. This leads us to a 2D velocity distribu- 
tion. We discuss the impact of the integration constant V', 
the disk radius fldisk, the normalisation factor A of vj/^^pec 
and the redshift factor a on the wind solutions. Table |^ 
summarizes the values of all other parameters which are 
kept constant throughout this paper. 

Table 1. Model parameters 



Model parameter 



Stellar radius 7?, 
Light cylinder Rl 
Injection radius rinj 

Core radius Rc 

Jet radius a of Eq. (120) 



Peak stellar magnetic field B, 

Dipole moment fj,. 

Stellar mass M* 

Stellar density n. 

Stellar Alfven velocity va,* 

Magnetisation a. 



Value 



3-Rq 

2250i?, 

1R» 

8R* 

1430i?. 

1000 G 

RlB, 

2Mq 

10^^ cm-^* 
690 km s" 
10"* 



P 7 



The distribution of the magnetic flux surfaces ^'(i?, z) 
and the corresponding magnetic flux tube functions $(i?) 
for each surface are shown in Fig. |l|. The critical surface 
^crit is the first open flux surface and reaches the asymp- 
totic jet radius i?jGt- There is no jet at radii larger than 
that given by ^crit- We choose our parameter a of Eq. 
(120) such that we get realistic jet radii of i?jet ~ (3— 4)i?i 
( [Fcndt fc Camenzind 199^ ). In our calculation the period 
of the star is P* — 1.14 days which is equivalent to a light 
cyhnder at Rl = 2250 R^. Because of the influence of 
the jet radius is shifted to bigger values than those given 
by a. We get i?jot = 7937i?, = 3.5i?L for a = U30R^. The 
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Figure 1. Top: Magnetic flux surfaces for Odisk = 2.5i?,, A = 10 and -0 = 0.1364. The contour levels are 
* = 1.01,1.00,0.96,0.93,0.86,0.78,0.75,0.71. The thick curve in the left panel represents the last closed flux surface 
and marks the beginning of the disk at adisk- Bottom: The corresponding flux tube functions $ 



core radius Rc = 8i?*. As injection radius for the plasma 
we choose rinj = IR*. 

In contrast to calculations done by Fendt & Camen- 
zind (1996) our analytical model is not limited by nu- 
merics, hence we can take a realistic weak magnetisation 
parameter of a* = 10~* which corresponds to observa- 
tions. Throughout the jet ct» is kept constant which is 
equivalent to a constant mass flow rate on each surface 



We normalize the radius R to the light cylinder Rl and 
the velocities to the Alfven velocity on the stellar surface 
VA.*, where v\ ^ = i?J/47rmpn,, and the flux surfaces to 
\I'crit- The density on the stellar surface is n* = lO^'^cm"^, 
corresponding to an outflow rate of Mjet — 10~^ M0yr~^. 

The flux tube function is a measure for the opening an- 
gle of the flux surfaces. Small $ correspond to large open- 
ings, constant values to a conical structure of \1/ ($ = const 
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represents a monopolar type flux surface). For small radii 
we can reproduce a stellar dipole which is mainly respon- 
sible for the strong decrease of the flux tube functions. 
The more the surfaces show a dipolar structure instead 
of a conical one, the longer is the extent of the region 
of decrease of the flux tube function. This is equivalent 
to an acceleration of the plasma. There is a deceleration 
of the plasma for the outer flux surfaces for radii larger 
than the disk radius. The corresponding ^-values then 
even catch up with those obtained by the more inward ly- 
ing surfaces. Nevertheless, their total plasma acceleration 
always exceeds that of the latter ones. The deceleration 
process coincides with a sharp bending of these flux sur- 
faces at the inner edge of the disk which is not found for 
the inner flux surfaces. From this we draw the conclusion 
that the inner dipole part of the star in our model is cru- 
cial for the plasma acceleration. The closer the flux surface 
is to the z-axis, the less the surfaces look dipole-like and 
show a conical structure instead. For intermediate radii 
the monopolar- type part of ^'^p'^'^ controls all flux surfaces, 
consequently all flux tube functions become constant. For 
large radii all surfaces coUimate. 



6.1. The cold wind 

A physical cold wind solution has to start at zero veloc- 
ity at the injection radius and, accelerating monotonically, 
pass through the Alfven point and the fast magnetosonic 
point. Figure ^ shows the wind solution along the flux sur- 
face = 0.89. The maximum poloidal velocity is reached 
when <I> is lowest. As a general feature it can be seen that 
the acceleration and deceleration of the plasma along dif- 
ferent flux surfaces varies. The acceleration for the outer 
surfaces is more efficient than for the inner ones and the 
asymptotic poloidal velocity u'^ is larger. Figure § shows 
the asymptotic velocity distribution over the cross-section 
of the jet. The velocity decreases for decreasing flux sur- 
faces. In our calculations, the mass flow rate per flux sur- 
face is conserved, ri{^) — const. 

The positions of the Alfven point xa and the fast mag- 
netosonic point xfm are plotted for several flux surfaces 
in Fig. ^. The resulting value for e is of the order of 
10~^. The Alfven point moves to smaller radii for de- 
creasing flux surface values This is equivalent to a 
smaller amount of angular momentum being transferred 
to the plasma. As expected for weakly magnetized flows 



(Fcndt & Camenzind 1996), the Alfven point is very close 
to the injection radius and far from the light cylinder. 
For successively decreasing surfaces the fast magnetosonic 
points move linearly inward towards smaller radii. For the 
innermost flux surfaces no determinations of xpM and 
i^crit are possible any longer. This means that for these 
surfaces the wind remains submagnetosonic. Wind solu- 
tions for flux surfaces lower than 50% of ^'crit start oscil- 
lating. 




R/R, 

Figure 2. Poloidal wind velocity Up along the flux surface 
* = 0.89 and for Odisk = 2.5i?,, V = 0.1364. The physi- 
cal branch starting with nearly zero velocity and passing 
through the Alfven and fast magnetosonic points is shown 
together with the unphysical branch which is retained in 
the subsequent flgures 



Concerning the question of inner wind solutions two 
different types of solutions are possible. Either we insist 
on every solution to be a critical one regardless of whether 
it is global or not, or we argue that the crucial feature of 
a wind is to be global regardless of being supermagne- 
tosonic. 

In the flrst case a flux surface can be found for which 
the acceleration is not sufficient for the plasma to reach 
the asymptotic regime (for 5" — 0.42 ^Pcrit, with the pa- 
rameter set fldisk — 2.5ii!,, A — IQ, ij} — 0.1364). In this 
case, the wind still reaches the fast magnetosonic point 
but is then so much decelerated that it closes very quickly 
afterwards (upper panel in Fig. ^, where m™'^ = u^pM- 
This implies that the inner part of the asymptotic jet is 
empty despite having a constant mass distribution for each 
^ at Xinj. Even though the flux surface extends to some 
asymptotic radius, the plasma flow stops somewhere on 
its way along 5*. Strictly speaking, the supermagnetosonic 
non-global solutions cannot be considered to be physi- 
cal solutions. This behaviour is due to the fact that we 
do not take pressure into consideration and have no self- 
consistent solutions. Especially close to the z-axis a wind 
can no longer be centrifugally driven and Lorentz-forces 
become negligible, hence pressure becomes more impor- 
tant. A mainly Poynting-flux driven cold wind (see Eq. 
(p3|)), where RB^ on the z-axis, will obviously not 
have any mechanism by which it can be accelerated away 
from the star. Only those winds can reach the asymptotic 



14 



E. Broitmoser & M. Camenzind: Collimated outflows of young steUar objects 




01 



R/R, 

Figures. Asymptotic velocities for the flux sur- Figure 5. Poloidal wind velocity Up. The solutions 
faces = 1,0. 96, 0.93, 0.89, 0.86, 0.82 and 0.78 and fldisk = are along the flux surface = 0.42 with Odisk = 
2.5i?*, ip = 0.1364 2.5i?*, V = 0.1364. The critical non-global solution for 

-E-crit = 1.0000005203 is shown in the upper panel, the sub- 
magnetosonic global one for i?ciit = 1.0000005215 in the 
lower panel. Here the physical solution is the lower branch 



On the other hand, we can argue that because of the 
small influence of centrifugal and Lorentz-forces near the 
z-axis a wind becomes less accelerated so that it can no 
longer reach the magnetosonic point. This wind flows only 
with submagnetosonic velocities, but reaches the asymp- 
totic jet radius. Hence we look for the smallest energy 
possible for which the wind can be global (Fig. ||). In 
this case the physical wind solution is the lower branch. 
We claim that the energy continues to decrease smoothly 
from one flux surface to the next as is the case for the 
critical solutions. This feature is found for the supermag- 
netosonic, global solutions and should be maintained. For 
small fj* (Fendt & Camenzind 1996) the principal accel- 




H/R, 

Figure 4. Magnetic flux surfaces 5* for Odisk = 2.5i?* and 
ip = 0.1364. The dots represent the position of the Alfven 
point, the crosses that of the fast magnetosonic point. The 
contour levels are "if — 1.0 to 0.64 in steps of 0.036 



eration of the plasma occurs at radii smaller than the 
Alfven-radius. There the plasma and the magnetic field 
corotate. This means that the acceleration is mainly due 
to the centrifugal force which decreases towards the z-axis. 
Consequently, less energy can be expected to be available 
near the z-axis to drive the wind. The final conclusion is 
that we have to investigate a hot wind in the future for 
describing the physics more consistently. 



radius for which the flux tube functions ^{R) decrease 
over a sufficiently long range. All other winds are not suf- 
ficiently accelerated. 



6.2. The influence of the disk radius ajisk 



We investigate the features of the wind for a constant 
value of the unnormalized critical surface ^f™™'''" (Ex- 
amples I-HI in Table 0; realistic disk radii are at about 
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Odisk ~ (2 — 3)i?*). This means we keep V'/^disk const, 
varying and Odisk at the same time, and impHes that the 
amount of flux passing through the hole between the star 
and disk and the asymptotic jet radius i?jct = 7937i?* re- 
main unchanged as well. Figure ^ shows the influence for 
three different choices of /psi — Odisk on ^f, $ and Up. For 
all three cases investigated, $ has the same value at Xinj . 
But the lower ip and fldisk are, the lower are the values of $ 
(lower ip corresponding to less steep flux surfaces ^), and 
these are reached at successively smaller radii (for succes- 
sively smaller adisk)- For all three parameter sets, the so- 
lutions for $ merge at a certain value of x, and thereafter 
track each other, remaining constant over the same length 
before they ultimately decrease to the same final value at 
X — Xjct- This behaviour results in the wind reaching the 
fast magnetosonic velocity at smaller radii for lower val- 
ues of ip and adisk- In addition, less energy is required to 
drive the wind, but with a smaller asymptotic velocity. We 
can conclude that a decrease of $ over a long distance is 
more important for obtaining large asymptotic velocities 
than a larger difference between $(i?inj) and the minimum 
$min (Table ^, Fig. ^. The step between obtaining a crit- 
ical, global solution or only a non-global one takes place 
between ^! = 0.42^'crit and ^ = 0.49^'crit- As already 
mentioned the wind solutions start to oscillate for flux 
surfaces which are lower than 50%^crit- This is true for 
all three parameter sets. Hence the exact determination 
of where this transition takes place is not always possible. 
As can be expected from an investigation of homogeneous 
solutions of the Grad-Shafranov equation, the parameter 
choices influence the physics most near the central source. 



Table 2. Three different parameter sets I, II and III for 
- 1 refers to ^fcrit = 1.00 





I 


II 


III 


"disk (-R.) 

(Scrit - 1) (10"") 

u^(km s-1) 


0.1364 
2.5 
1.232 

444 


0.2727 
5.0 
1.234 
445 


0.6000 

11.0 

1.413 

482 



6.3. The weighting between 5'''°'° and 

Now we investigate the influence of the constant A of Eq. 



(120) on the wind solutions. All other parameters are kept 
fixed and all solutions are for ^ = Vl/crit- Figure |^ shows 
the characteristics of the flux surfaces and the flux tube 
functions for v4 = 5, 10, 120 and 240. A larger A-value has 
the same effect as increasing the integration constant -0: 
more flux can escape to infinity and less is forced to pass 
through the hole between the star and the disk. The po- 
sition of the asymptotic radius is shifted to higher values 
for decreasing A. For example for A = IQ the critical flux 



1* 




R/R, 




Figure 6. Flux surface (upper panel) ^' — 0.86 and flux 
tube function (middle panel) for adisk = H, = 0.6 (dot- 
ted curve), adisk = 5, '0 = 0.2727 (thick curve) and for 
adisk = 2.5, = 0.1364 (dashed curve). The lower panel 
shows the poloidal velocity for adisk = H, = 0.6 (thick 
curve) and adisk — 2.5, psi — 0.1364 (thin curve) 

surface will collimate at radii about flve times larger than 
for A = 120. This is due to the larger impact of on 
the total solution which will shift the asymptotic radius 
to a higher value. The higher the value of A, the greater 
^(ccinj) is and its rate of decrease for small radii, and the 
closer to adisk does the critical surface reach its minimum 
$-value. This results in a wind reaching higher velocities 
near the star, and the fast magnetosonic velocity becomes 
larger for larg er A - values (Fig. In contrary to the re- 
sults in Sect. 3.2, the maximum velocity is higher when 
the flux tube function decreases more rapidly but for a 
shorter distance (compare A — 240 to A= 120). Here the 
stronger decrease of $ for A = 240 overcomes the fact that 
it only decreases for a shorter distance. After $ reaches 
its minimum value, the flux tube function increases for all 
A resulting in a deceleration of the wind. But for inter- 
mediate radii the features of the flux tube function vary 
depending on A. For small A, <I> becomes constant and 
even decreases in the collimation region. For larger A the 
flux tube function is increasing. Again there is a decrease 
in the collimation region. Nevertheless this continuing in- 
crease of $ for intermediate radii is crucial for the wind 
to be able to reach the asymptotic radius. For A — 120 
the increase in the flux tube function is not so steep that 
a global wind solution could no longer be found. On the 
other hand, for A — 240 no global solution any longer ex- 
ists despite the flux tube function having qualitatively the 
same appearance as that for A = 120. 
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H/R, 

Figure 7. The flux surfaces (upper panel) and the flux 
tube function (lower panel) for A — 240 (long-dashed), 
A = 120 (sohd), A ^ 10 (dotted) and A = 5 (short- 
dashed). All calculations are done for the critical flux sur- 
face, 4* = 1.00, Odisk — IIR* and ip = 0.6 



6.4- The role of the flux integration constant tp 

The integration constant tp represents the ability of the 
homogeneous solution for the flux surfaces to open and 
to let flux escape to infinity. The higher ip is, the higher 
the bending of the dipole-like shaped field line after its 
opening near Cdisk- For the infinitely small disk this occurs 
nearly at ^ = 0. For -0 = Eq. (116) will not contribute 
to radii R > Odisk and the magnetosphere is forced to go 
through the hole between star and disk. We get a dipole 
which is squeezed by the inn er ed ge of the disk. But due 
to the special solution (Eq. ( |120D ) plasma can escape to 
larger radii. Then the shape of the outer flux tube function 
is only determined by 5'''?'='=, consequently i?jot = a. For 
asymptotic radii the flux tube function increases up to the 
final jet radius, which is in contrast to the solutions we get 
for tp = 0.6 (see the result for A = 10 in Sect. 6.3), where 
the flux tube functions decrease for asymptotic radii. For 
both cases the wind can reach the asymptotic region. 

The influence of ^ on the shape of the magnetic flux 
surfaces is discussed now for five different values, as listed 
in Table || and plotted in Fig. ^, and for the case of the 
critical surface ^'crit- For ip — 2.0 the flux surface opens 
at Rq « 4.5i?* — fldisk and then keeps very close and 
flat along the surface of the disk. For tp — 2.5 the flux 
surface dips down to z = at i?o ~ 4.0i?* < Odisk, then 
rises instantly, but only to a height much smaller than 
that which the dipole reaches as a maximum, and becomes 
constant about li?* later. Increasing ■0 to 3.5, Rq is further 



R/R, 



Figure 8. The thin curve shows the poloidal wind for A = 
120, the thick curve that for A — 240. It is 4'crit, Odisk = 
Hi?, and ^ = 0.6 

Table 3. Five different values for the integration constant 
C and the related radius Rq, at which the critical flux 
surface opens 





I 


1—1 


III 


IV 


V 




2.0 


2.5 


3.5 


4.5 


10 


.T.unnorm 
^crit 


1.32 


1.37 


1.43 


1.48 


1.64 


Ro (R,) ~ 


4.5 


4.0 


3.5 


3.0 


< 3.0 



lowered to about 3.5i?,. Afterwards the flux surfaces rise 
very quickly without having a constant part parallel to the 
X-axis. These characteristics occur and are enhanced for 
ip — 4.5 and tp = 10. For these examples the normalisation 
factor A = 10 is used. For smaller values of A the influence 
of the homogeneous solution is so great that the flux 

surfaces would collimate at much smaller radii (e.g. for 
A^ 10 they coUimate at about 8000 i?,, for A = 120 at 
about 1650 R^). 



6. 5. The role of gravity 

We next consider the role of gravity on the wind. Note 
that in principle the vacuum solution \I>hom should now 
be derived from the Grad-Shafranov equation (Eq. ([78[)) 
including the gravitational redshift a. But it can be shown 
that the influence of a for the solution is negligible (see 
Appendix C). 

The importance of the redshift factor a on the wind is 
seen at small radii. A wind can only emanate from the star 
if (£■ — 1) > 0. From Eq. (|5^) we recognize that for a < 1 
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10 



Figure 9. The critical flux surfaces for = 2.0, 2.5, 3.5, 
4.5 and 10 and Odisk — 4.5i?*, A = 10. All plots have the 
same scaling and are shown from i? = (0 — 6)B„ 



solutions with positive energy exist only for Up{xiaj) > 
which cannot be satisfied by a cold wind (p* ^ e). This 
also means that the critical energy can no longer be calcu- 
lated by requiring mQ^Xinj) = (see Eq. ([z^)) because this 
implies Wp(a:inj) = 0. We can conclude that the cold hmit 
is not a good approximation for our problem. The devia- 
tion from a = 1 by taking gravitation into consideration 
is of the order 10^^. That is the same order of magnitude 
as for the expression e/(l — e). 

Figures ^ and [ill demonstrate the main differences 
between a wind calculated for the non-gravity limit {a — 
1) and the case including gravity. Very close to the star 
(a < 1) the non-Newtonian wind cannot start with zero 
poloidal velocity. We claim that the ratio of sound speed 
to Alfven velocity on the stellar surface is Cs^^,/va.* = 
Up{xinj)/vA,* < 1. A sound speed of Cg.* — 100 km s^^ 
in a stellar corona can be assumed. As can be seen in 
Fig. |Tl| this ratio is fairly higher in our calculation (about 
0.73 for e = 1.26 10~^). Instead of determining the critical 
energy and e for a given injection radius x-mj , we can choose 
a higher value for e and still obtain a critical solution 
starting at about the same injection radius but with a 
smaller velocity. Along = 1 reasonable results can be 
obtained up to e = 1.8 x 10~^ and Cs.^,/va.* = 0.55. This 
value is still too high. 




H/R, 

Figure 10. Magnetic flux surfaces for Odisk = Hi?* and 
^l> = 0.6. The inner dots represent the positions of the 
Alfven points and the outer ones those of the fast mag- 
netosonic points, for the wind without gravity. The inner 
crosses show the positions of the Alfven points and the 
outer ones those of the fast magnetosonic points, for the 
wind with gravity. The contour levels are ^E* — 1.0 to 0.64 
in steps of 0.036 and * = 0.49, 0.57 




R/R, 

Figure 11. Poloidal wind velocity Up. The solution is 
along the flux surface = 1. The thin curve is calcu- 
lated with a — 1, the thick one with a < 1. The inner 
disk radius is Odisk = Hi?* and V' = 0.6 
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The wind reaches the asymptotic radius but for lower 
velocities due to the deceleration caused by gravity. As a 
consequence, the value of the critical energy is smaller. In 
comparison with the case of a constant a, the Alfven-point 
for a radial dependent a{r) for all flux tube functions ^ is 
reached at smaller radii, and the fast magnetosonic point 
at larger radii. Their shapes are nonetheless qualitatively 
similar. For a calculation with or without gravitation, the 
Alfven-points as well as the fast magnetosonic points move 
linearly inward with decreasing Figure |lo| also shows 
that for a < 1 a critical wind solution cannot be found for 
flux tube functions as small as for a — 1, i.e. the Alfven- 
point and the fast magnetosonic point can no longer be 
determined. The position of the Alfven- surface and the 
fast magnetosonic surface stop rising linearly on the sur- 
faces for which "if ~ 0.57^crit- They appear to run parallel 
to the flux surfaces, i.e. the z--axis, a result also obtained 
by Fendt & Camenzind (1996). For both calculations the 
fast magnetosonic points are at radii smaller than the disk 
radius. For the parameter set (I) of Table |2| a wind calcu- 
lated along the flux tube function 'i' — 0.57 and including 
the gravitation of the star can no longer reach the asymp- 
totic region. 

6. 6. Derived wind parameters 

Figure ^ shows the poloidal velocity for '5 ~ 0.89 and 
the corresponding Alfven-Mach number M, the poloidal 
current /, the toroidal velocity and the density n. A 
constant poloidal current / is equivalent to a force-free 
magnetosphere. But the current strength decreases to- 
wards larger radii by a factor of about 7 — 8. The rela- 
tion between Up and is a subject of interest when the 
observed emission lines are to be interpreted. There is a 
high velocity dispersion if Up and are of about the same 
order of magnitude, hence the spectral lines are strongly 
Doppler broadened. In our model these two velocities are 
of about the same order of magnitude for small radii. At 
the Alfven point, one finds Up/u^ ~ 0.7, at the fast mag- 
netosonic point Up/ucf, ~ 3.7 and at the asymptotic radius 
Up/uif, ^ 2.2 X 10^. The plasma density between the stel- 
lar position and the asymptotic jet radius decreases by 
several orders of magnitude. 

7. Conclusions 

In this work we presented an analytical model for the mag- 
netic flux surfaces which reproduces the main features of a 
protostellar wind that coUimates at large radii. The com- 
plete expression for the flux surface consists of three differ- 
ent parts which determine the solution for small, interme- 
diate and large radii. For the region near the star a stellar 
dipole field is assumed and calculated by extending the 
computation of Kundt & Robnik (1980). We assume the 
existence of an infinitely thin disk perpendicular to the ro- 
tation axis at a distance of a few stellar radii (odisk = Hi?* 



for most examples). The stellar dipole field turns out to 
be essential for the acceleration of the wind to supermag- 
netosonic velocities, which occurs very near the star. For 
intermediate radii the magnetosphere is nearly spherical 
and may be described as monopole. CoUimation of the 
field fines is achieved at radn i?jot « (3 — 4)7?^. This re- 
gion is characterized by magnetospheres similar to those 
for force- free solutions (Appl & Camenzind 1993b). 

We investigate the dynamics of a wind in these colli- 
mated magnetospheres for axisymmetric and stationary 
flows, obtaining a 2D picture of the wind. As a first 
approximation we consider a cold wind, neglecting the 
plasma pressure. The physical solution must start at zero 
poloidal velocity and pass smoothly through the Alfven 
point and the fast magnetosonic point. The resulting wind 
solutions reach the jet radius in the Newtonian limit. The 
asymptotic jet velocities of about a few hundred km/s 
are comparable to those observed and depend mainly 
on the choice of the magnetization parameter cr*. The 
asymptotic velocity of the outflowing plasma decreases 
the closer it is to the jet axis. Realistic asymptotic jet 
velocities are obtained for the parameters cr* — 10~^, 
a = 1430i?,, A = 10, = 0.6, Odisk = H-R* and = 482 
km s^^ for ^'crit. The influence of the various parameters 
of the model on the wind solutions is discussed in more 
detail in the text. If we include the effects of gravity by 
means of the redshift factor which accounts for the gravi- 
tational force of the central star, solutions must start with 
a non- vanishing velocity. This is equivalent to a rise in the 
sound speed up to about cs.*/va.* < 1- Such a high ratio 
of sound to stellar Alfven velocity could only be explained 
by temperatures near 10^ K. By contrast, shock heating 
due to plasma accretion onto the star allows only temper- 
atures up to 10^ K. Hence an additional source of heating 
is required. Another explanation might be that a hot wind 
is essential to obtain a solution. 

For plasma flowing along magnetic flux surfaces which 
are close to the z-axis, critical wind solutions can no longer 
be determined. Either there is no global solution branch 
extending to the asymptotic radius of the flux surface, or 
there is a global plasma flow which remains submagne- 
tosonic throughout the entire flow. The last two points 
clearly demonstrate the need to calculate a hot wind in- 
stead. 

The velocity field throughout the jet and the related 
density and temperature obtained here can be used for 
future work as input to calculate forbidden emission line 
profiles of T Tauri Stars ( ^olf fc Bohm 19931 ). These for- 
bidden emission lines are a powerful means of probing the 
outflow of Classical T Tauri Stars and of studying the link 
between the ejection and accretion processes (Breitmoser 
& Camenzind 2000). 
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Figure 12. Poloidal wind velocity Up, Alfven Mach number M, current J, toroidal wind velocity and density n. 
The solution is along the flux surface ^ = 0.89. All curves are calculated for a = 1, adisk = Hii* and ip = 0.6 
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Appendix A: Homogeneous solution of the 
Grad-Shafranov equation 

A homogeneous solution of the Grad-Shafranov equation, 
written in cyhndrical coordinates. 



(Al) 



d^A^ d^A^ IdA^ _ 
dR^ dz^ R OR i?2 

can be obtained by separation of variables 

/>OC 

A^{R,z)= I b{k)Ji{kR) exp{-k\z\)dk (A.2) 
Jq 

given by suitable boundary conditions. The magnetic flux 
surfaces are given by 5" = RA^p, and Ji{kR) is a Bessel 
function of first kind and first order. We make the follow- 
ing assumptions: 

— The plasma is perfectly conducting, consequently no 
field lines can penetrate the disk and the disk is 
screened by induced surface currents. 

— The disk is infinitely small and exists therefore only at 
z = 0. 

— The dipole is perpendicular to the disk. 



The problem is axisymmetric hence the in- 
duced vector potential A is given by A^p alone 



( |Kundt fc Robnik 1980| ). 
The two boundary conditions are 

- Br{R, z = 0) = for i?, < i? < Odisk 

- Bz{R,z = 0) = for R > adisk- 

For the most general case of open flux surfaces the 
second condition, in contrast to Kundt & Robnik (1980), 
reads 



<iP°l^(i?,0) 



Const 
R 



b{k)Ji{kR) dk. 



for R > fldisk- 



(A.3) 



In addition to their result we get a part proportional to 
the integration constant Const 



b{k) = 



t''^,h/2{kt)dt 



d 



Const 



= Const 



s\n{kt) cos(fci) ' 



-.dT 



kt^ 



t 



dt 



(A.4) 



and 

- 



Const -J-t^-2i\z\t + R'^ 

Im — 

R t 

Const / 1 

■ 7 



(A.5) 



sm f \p 

For our further calculation we use '0 = Const adisk7''/2/i*. 

Appendix B: Magnetic structure of collimated 
magnetospheres 

For the magnetic fields B = 5^°'" + B"P'='= we get 
1 



Br — 



Bb 



sin do ' 
1 



B 



r sin 6 dr ' 

hom _ cos 



arctan(7) -f - (1 - 5 + C- 



(B.l) 
(B.2) 

(B.3) 



^hom ^ — arctan(7) + ■ ^ ' ' 



sin^e 



+C^(l-5(l+7'))), 
sm 



+ 2 cos 6 1 



where 



r^i?^ arg 
2 sin 61 1 
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Appendix C: The Grad-Shafranov equation 

including the gravitational redshift 



In comparison to Eq. (A.l) the Grad-Shafranov equation 
including the gravitational redshift (Eq. d78)) reads 



with 



R 



B.dA 1 
R^' dR ^ i?2 



B = 



(C.2) 



Eq. (C.l) differs from Eq. (A.l) in the dimensionless ex- 
pression B/R. This has a maximum value at R = R^. An 
estimate shows that B/R^ « 10~^^ is negligible for the 
calculation of 'J. 



